The Ratio of Gluon Distributions in Sn and C* 

T. Gousset and H.J. Pirner 
Institut fur Theoretische Physik, Universitat Heidelberg 
Philosophenweg 19, D-69120 Heidelberg, Germany 



Abstract 

We calculate the ratio of gluon densities, G Sn (x) / G c (x) , for 0.01 < x < 0.1, 
from the new high statistics data on F^/F^ taken by the NM Collaboration. For 
small x, the shadowing in the gluon distribution is about equal to the shadowing 
of quark distribution. The antishadowing in the gluon distribution, however, is 
roughly 10%. We also compare with the ratio difference R = <jl/&t horn Sn and 
C. 
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The investigation of nuclei with deep inelastic lepton scattering has produced very 
accurate data on nuclear structure functions F 2 A (x). Especially the measurement of the 
ratio of tin to carbon structure functions has been carried out recently with very high 
statistics ]]J. It shows shadowing for x < 0.05. Antishadowing (< 2 %) is visible around 
x = 0.1, i.e. at distance scales of r = 1/xMjy rs 2 fm corresponding to the average distance 
between nucleons. In the interval 0.15 < x < 0.8 the nuclear valence quark density is 
reduced. By measuring three different yU-beam energies (120, 200 and 280 GeV) NMC has 
also been able to determine the weak In Q 2 dependence of the structure function ratio. 
The aim of this note is to use the leading order evolution equations to estimate the little 
known underlying gluon densities. 

There are older data from inelastic J/ij) production at z < 0.85, p\ > 0.4 (GeV/c) 2 , 
which determine G Sn (x)/G c (x) ~ 1.13 ± 0.08 in the interval 0.05 < x < 0.15. Hadronic 
production of J/if>, ip' and T in pA or it A collisions [^] has been analysed with the 
gluon-gluon fusion process to get G A (x 2 )/G p (x 2 ) in the region where the nuclear gluon 
momentum fraction x 2 is small [f|, [§]• In the small x 2 region, however, the second gluon 
is as soft as the third gluon which is necessary to form the J/ip 3 5"i-state. Therefore one 
expects final state interaction effects in the small x 2 region [El. 

Currently there is a discussion to use nuclear beams in HERA. The study of the 
gluon distribution in nuclei at small x would be one of the outstanding new opportunities 
to investigate nuclei on the parton level. Also the forthcoming heavy ion experiments at 
RHIC and LHC need a good understanding of the nuclear parton distributions to calculate 
the cross sections for hard processes initiated by nuclear collisions. 

Measurements of the gluon distribution in nuclei give experimental windows viewing 
the partonic structure of nuclear binding. Very little is known about the role gauge fields 
play in the nucleus. This has been studied for the case of abelian QED in a model 
where the hydrogene atom replaces the nucleon. In this model, one computes the photon 
density of the hydrogene molecule, G H2 (x), and compares it to the free atom photon 
density, G R (x). The electron orbits of the hydrogene atoms in the molecule are polarized 
and modified by the electron exchange interaction leading to a suppression of photons 
at small x. At the momentum corresponding to the relative distance of the two protons 
a small antishadowing peak is visible |7|] . In analogy, gluon antishadowing may indicate 
the relative distance (xMn)^ 1 ~ 2fm between the centers of nucleons, which act as a 
source for color fields. A covalent binding of quarks may manifest itself as a density 
dependent lack of long range gluons at x < 0.1 similarly to the deformation of the photon 
cloud in the hydrogene molecule. In addition, in non-abelian QCD one expects at small 
x processes where gluons from different nucleons overlap and merge || . Both effects have 
also an interpretation in the nuclear rest frame in terms of absorption of various hadronic 
components of the photon in the nucleus. 

The strong theoretical interest, the excellent new data of NMC and the starting discus- 
sion of HERA with nuclei motivated us to extract from the available data the maximum 
information on the gluon distributions in nuclei. 



2 



As usual, we define the Bjorken variable relative to nucleon kinematics 
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with Mtv the nucleon mass, v the virtual photon energy in the laboratory and Q 2 = —q 2 
its virtuality. The Leading-Log (LL) approximation reads (for a review see Ref. 0) 

F 2 (x,Q 2 ) = x^2e 2 qi(x,Q 2 ), 
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where qi and G are the quark and gluon densities in the nucleon or in the nucleus (in 
the nuclear case, we shall consider distributions per nucleon). The splitting functions are 
denoted by P qq and P q c- The sum runs over active flavors of quarks and antiquarks and 
ti is the corresponding charge. It is understood that distributions, structure functions, as 
well as as are Q 2 dependent. As we will see, for x sufficiently small, the second term in 
the evolution equation of F 2 is small in comparison with the first. Furthermore, the first 
term can be approximated in a simple way to easily extract the gluon distribution from 
the evolution of F 2 flQfl . In the nucleon case, the LL approximation becomes insufficient 
at very small x. Next-to-LL (NLL) corrections are about 10 — 20% at x = 0.01. For our 
purpose of extracting ratios in the range 10~ 2 — 10 _1 , however, the LL approximation is 
accurate enough. In the nuclear case, additionnal merging terms are expected to play an 
important role in the evolution equation with the increase of the gluon density. These 
non-linear effects are expected to appear at very small x (< 10~ 2 ) ||. 

Let us consider the first term in Eq. (0). It can be written in detail as 
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We can replace the lower bound of the integral, x/A, by x/A — > 0, if in addition we 
impose G(u) = in the unphysical region, u > A. Because the gluon density is completly 
suppressed for x > 0.5 this replacement is an accurate approximation in the small x region 
for any reasonable ansatz for G(x). We then expand 
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around y — 1/2, 
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Carrying out the integration, odd terms in the expansion do not contribute and one gets 
f Q dy(y 2 + (1 - yf)g{x/y) = ^ g(2x) + ^ [(2x) 2 g"(2x) + 2(2x)g'(2x)} +■■■ (4) 
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Figure 1: The correction —A(x) for different Q 2 in the nucleon case. The NMC fit to the gluon 
density |11], evaluated at 2x, is also plotted for comparison (the Q 2 for this density is 7GeV 2 ). 



As found in [jnj, this expansion is efficient for standard gluon distributions and one can 
approximate 

f 1 dy(y 2 + (1 - yf)g{x/y) « \ g(2x). (5) 

JO 6 

One can effectively check that for g{u) = K u~ A , the right hand side of Eq. (^) becomes 

!*<w-(i + ^ + -). 

showing that Eq. (|5|) is accurate for < A < 1, with an error smaller than 5%. 
We now consider the second term in Eq. (0) 

A(,, Q») = x (p„ 3) = 1 /; * C 1 ^ 2 '^'- 2 ^' 4(1- B )F,(,/ V) ) ; 

V x J 6 Jx/a \ 1 — y 2 J 

(6) 

One can compute A(x, Q 2 ) from F2 data. We show in Fig. [I] the behavior of A as a 
function of x and Q 2 in the case of the nucleon. The behavior of the gluon distribution 
is also plotted for comparison. For x ~ 0.2, A(x) is of the same order of magnitude as 
2xG(2x) and indeed cancels the gluon contribution in Eq. (0) leading to the true scaling 
of F 2 . Due to the steep rise of G, A becomes a small correction as soon as x < 0.1, thus 
implying a simple relation between F 2 evolution and gluon density. 

In the following, we shall take into account both contributions to the evolution equa- 
tion, Eq. (H), i.e. 
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We can now turn to the study of the new NMC data on tin to carbon ratio of structure 
functions. Namely, measurements of 
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have been performed JFJ. To be more precise, let us mention that the ratio, f\, has been 
measured for fixed x over a certain range of Q 2 , so that the reported value is an average 
of ratios corresponding to an averaged (Q 2 ), where (Q 2 ) varies from one x to the next. 
The evolution of ratio with InQ 2 has then been computed, giving the "slope", /2, at x 
and (Q 2 ). 

Using Eq. (ffl), it is easy to relate the tin to carbon ratio of gluon distributions 
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In this analysis, we use the fact that the ratio of carbon to deuteron structure functions 
shows practically no Q 2 dependence |l2 |, so that we can drop the second term in the right 
hand side of Eq. (|TTD, and use for f 3 the deuteron data of [ pT[| . 

The extraction corresponding to Eq. (|ToD is given in Table |l] for each x together with 
the corresponding (Q 2 ). We plot the ratio r(x) in Fig. |2| together with fi(x). One 
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Table 1: The ratio of tin to carbon gluon density for different x and corresponding {Q 2 ). 5r is 
the uncertainty on r coming from the statistical error on the slope /2- 
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Figure 2: The ratio r(x) = G Sn (x)/G c (x) of tin to carbon gluon density as a function of x 
together with the ratio of structure function, fi(x) = i^ n (a;)/ F^ix). The statistical error on 
fx is less than 1% in the whole range of x. The box represents the extraction of r from J /if) 
electroproduction data M (see text). 

observes a signal for a stronger antishadowing in the gluon case (~8%) than in the quark 
case (~1%). In the shadowing region, there is no evidence for a difference between quarks 
and gluons. It is, however, presently impossible to draw a firm conclusion in this region. 
It turns out that the uncertainty in the right hand side of Eq. ([TDD comes predominantly 



from the error in f% so that, in Table [I], we report the uncertainty 
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where A/ 2 is the experimental statistical error on f 2 . The systematic error is in the 
process of being assessed. 

There are also theoretical approximations. Neglect of NLL corrections together with 
the use of Eq. (|5|) lead to an approximation presumably better than 20% for an absolute 
extraction of gluon density. Because the ratio r(x) remains close to 1 in the whole range 
of x being considered, they are hopefully even far better for the ratio we are evaluating. 

The ratio r has also been extracted from J ftp electroproduction || in the range 0.05 < 
x < 0.15, with the averaged result (r) = 1.13 ± 0.08. In Fig. 0, we show this value as a 
rectangle the sides of which are the x range of the data and the errorbar in (r). One can 
see that the two different extractions are in good agreement. 

Using the result on r, one can estimate the consequence on AR = R Sn — R c , the 
difference of the ratios of longitudinal to transverse virtual photon cross sections. In 
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principle measurement of ol at small x provides a direct information on G(x). This 
measurement is, however, difficult, therefore experimental analysis of ol is still going on 
in the small x range. Here we examine the problem by going in the opposite direction, 
that is to see whether the extracted ratio of gluon distributions, r(x), is compatible with 
the NMC measurement of AR IJ. 
From the expression of R, 

R = 2 ~ 1 — ; (12) 

Fi — Fl F2 — Fl 

assuming AR <C R (see below), one finds 

/ 1 pSn \ 

ARkR(1+R)(--± t -1). (13) 



At next to leading order, Fl is given by || 

F L (x) = ^x 2 jf* ^ (^F 2 (y) + 2 £ e?(l - x/y) yG(y)) . (14) 

At small x, the integrals can be approximated in a way similar to the method leading to 
Eq. (D. One has § 

Fl{x) « g ^F 2 (2x) + ^2.5*G(2.5*)) . (15) 

The iVterms are smaller than the G-terms, but they are not negligible in the whole x 
range. A detailed analysis would require a better knowledge of absolute values of F 2 's and 
G's, so we restrict ourselves to an estimate for AR by dropping the F 2 -terms to simply 
end with 

AR r(2.5x,Q 2 ) 

-1. (16) 
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Simple parametrizations of r and fi shows that this quantity varies very weakly with x in 
the range 0.01 < x < 0.1. It is about 0.08±0.01. R is about 0.1-0.3 for the nucleon case 
in this x-range (and for the Q 2 of the data), so that one verifies the above statement that 
AR -C R, and one gets 0.01 < AR < 0.03. This is compatible with the new data from 
NMC Q. We note, however, that the present precision of the data on AR unfortunatly 
does not give much constraint on the extracted ratio r(x). 

The present analysis can be compared with previous analysis ||, || of J/ ^-production 
in hadronic collisions 0. In an older attempt [[|, the experimental data have been 
all accounted for by gluon shadowing. The resulting ratio, G Sn /G c ~ 1 + (0.021 ± 
0.001) lnx2 In 118/12, would be much below our determined ratio in Fig. for x > 0.01. 
Recently, the same data have been reanalysed |J on the basis of a Fock state decomposi- 
tion of the charmonium states with an essential contribution from the \(cc)gg) component, 
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supposed to give the interaction of the precharmonium state soon after formation. This 
new analysis can explain the FNAL data without shadowing of the gluon distribution at 
variance with our analysis. 

The shape of r(x) obtained in this paper has been anticipated in [T3| on the basis of 



the fulfilment of the momentum sum-rule. In order to get a momentum fraction for the 
gluons ~2% higher in the Calcium nucleus than in a free nucleon, shadowing in the gluon 
distribution at x < 0.05 was predicted to be compensated by a noticeable antishadowing 



around x = 0.1. We can infer that the general feature discussed in |L3[ are in qualitative 
agreement with the data. The sum-rule fraction for the shadowing+antishadowing region, 
x < xo ~ 0.25, is 

Jo' 25 dx (r(x) - 1) xG(x) 
Jq dxxG(x) 

in tin to carbon. This value is close to the estimated enhancement of the gluon fraction 
in Ca to D. Let us note for completeness that the contribution from x > 0.25 to the 
momentum fraction carried by gluons is not totally negligible, in spite of the small gluon 
density, G(x), at large x. Typically, a mean depletion of 10% in the interval 0.25 < x < 0.5 
(EMC effect for gluons) diminishes 7^ from 2% to 1%. 

The range of x values accessible via the NMC ^-measurements 0.01 < x < 0.2 is 
sufficient to estimate the importance of shadowing for the minijet production rate at 
RHIC. One needs, however, much smaller x information, down to x ~ 10~ 4 , for heavy 
ion collision at LHC. To cover this range, HERA with nuclear beams would be very 
helpful. The scaling violation analysis shown here may have difficulties to be used since 
the standard evolution of the F2 structure function, even with NLL corrections included, 
may be modified by non-leading twist terms coming from the merging of partonsQ (EL It 
is therefore important to extend the inclusive measurements of F 2 to jet, open charm and 
charmonium productions. 

We gratefully acknowledge two members of NMC, Antje Brull and Andreas Mucklich, 
for providing us the data and related informations. We have also benefitted from discus- 
sions with participants of the working group on nuclei at HERA. 



1 Let us mention that this non- linear effect rather than being a limitation to gluon density analysis is 
to be seen as one of the very challenging effects to be studied at HERA. 
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